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For a system of non-interacting electrons scattered by neutral impurities, we derive a modified 
Boltzmann equation that includes quasiparticle and virial corrections. We start from quasiclassical 
transport equation for non-equilibrium Green's functions and apply limit of small scattering rates. 
Resulting transport equation for quasiparticles has gradient corrections to scattering integrals. These 
gradient corrections are rearranged into a form characteristic for virial corrections. 



I. INTRODUCTION 

In the first paper of this seriesEI (referred as Paper I) , we 
have discussed interplay of quasiparticle and virial correc- 
tions for scattering by neutral impurities with resonant 
levels. We have found that both corrections are of the 
same order and tend to mutually compensate. Accord- 
ingly, one should either include both kinds of corrections 
or neglect both of them. Separate quasiparticle or sepa- 
rate virial corrections lead to overestimates of impurity 
effect on basic physical quantities like dc conductivity or 
screening length. 

In Paper I we have used an intuitive modification of 
the Boltzmann equation (BE), Eq. (1-39) [Eq. (39) of 
Paper I] . Long time experience with quasiparticle correc- 
tions shows that intuitive approaches whatever convinc- 
ing are far from to be reliable as the wave-function renor- 
malization factor often emerges in an unexpected way. To 
become trustworthy, Eq. (1-39) has to be recovered from 
quantum statistics in a very systematic manner. This is 
the aim of this paper. 

The intuitive modification of the BE we would like to 
arrive to reads 

21 + z ^. d J- _ d A d J_ = -L 

dt m dr dr dk t 

+ 7¥/^ (|p| - |fc|)/(P ' r ' i - At) - (1) 

This is equation (1-39). It applies to system of non- 
interacting electrons scattered by point impurities with 
resonant levels. Here, the quasiparticle distribution / 
is a function of momentum k, coordinate r and time t. 
Lifetime t, collision delay A t and wave-function renor- 
malization z are functions of momentum, potential <f> de- 
pends on coordinate and time. The collision delay A t 
(given by the energy derivative of the phase shift) makes 
the scattering-in integral non-local in time. This non- 
locality represents virial corrections. The quasiparticle 
corrections are covered by the wave-function renormal- 
ization z. 



So far, quasiparticle and virial corrections have been 
studied only separately using different theoretical tools. 
Now we briefly review previous studies to identify which 
tool is better suited for unified theory. 

A. Virial corrections 

A need to derive virial corrections to a quantum trans- 
port equation has been felt for a while. The progress 
achieved ia. this direction can be represented by Snider's 
equationHO Using the BBGKY hierarchy, Snider has de- 
rived a quantum transport equation for the reduced den- 
sity matrix (Wigner's distribution). Snider's equation 
is sufficiently general to describe various corrections be=. 
yond the BE, however, as noticed by Laloe and Mullin,EJ 
all corrections enter the scattering integral, while quasi- 
particle corrections should appear rather as corrections 
to drift of single-particle excitations. 

Snider's equations includes another item alien to the 
quasiparticle picture: the reduced density matrix. When 
the quasiparticle corrections are in game, the transport 
equation should deal with the quasiparticle distribution 
not with the reduced density matrix. In the reduced den- 
sity matrix, contributions from free excitations (quasi- 
particles) and from correlations (off-shell motion during 
dressing processes) are mixed together. The drift of free 
excitations can be described by simple quasiclassical tra- 
jectories while correlations require quantum mechanical 
treatment. To be able to make efficient approximations, 
these two types of motion has to be separated. So far, 
the theory of transport based on the reduced density ma- 
trix is missing a tool that would furnish us with such a 
separation. 

B. Quasiparticle corrections 

The second group of papers is a large variety of stud- 
ies recovering the Boltzmann equation with quasiparti- 
cle corrections in the form visualized by Landau, e.g. 
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Refs. |-||. On the other hand, none of these studies 
touches virial corrections. One expects that the virial 
corrections should emerge from systematic quantum- 
statistical approaches to the Landau theory provided one 
does not loose them making unjustified approximations. 

^From the intuitive modification of the BE (0), we can 
identify three groups of approximations we have to avoid: 

• Most of authors limit their attention to the scat- 
tering described within the Born approximation (as 
we have also done in Ref. |lj). The corresponding 
T-matrix is just an impurity potential that does 
not cause any collision delay. Similarly, the self- 
consistent Puff-Whitefield approximation of the 
electron-phonon scattering (the Migdal approxima- 
tion with non-self-con&istent phonons) studied by 
Prange and Kadanoff,B gives no virial corrections 
since the interaction vertex has no internal elec- 
tron dynamics like the T-matrix within the Born 
approximation. 

• Prange and Sachs0 have studied the electron- 
electron and electron-phonon scattering within the 
fully self-consistent single-loop approximation. The 
screened Coulomb interaction is in general a com- 
plex function of energy so that it yields some kind 
of virial corrections. Prange and Sachs have used, 
however, the static approximation within which the 
screened Coulomb potential looses energy depen- 
dence and becomes real. In this way, the virial 
corrections are lost. 

• Danielewicz0 and Botermans and MalffietEl have 
used the two-particle T-matrix for nucleon-nucleon 
interaction that definitely includes virial correc- 
tions, however, they have neglected gradient con- 
tributions to the scattering integral. From the non- 
local character of the scattering integral in the intu- 
itive BE (0), one can see that the virial corrections 
are proportional to gradients of the quasiparticle 
distribution, see (||). Neglecting gradient contribu- 
tions one looses the virial corrections. 

In all the three cases, one can go beyond these approxi- 
mations. To recover Eq. (^), we have to use the T-matrix 
for the impurity scattering and keep the gradient correc- 
tions to scattering integrals. Similar treatment— has lead 
to virial corrections for two-particle scattering.E_3 

C. Gradient corrections to scattering integrals 

The gradient corrections to the scattering integral are 
of central importance in our treatment. In the intuitive 
BE ([!]) , the only gradient contributions to the scattering 
integral come from the virial corrections, 

d A + _d±dl = J_ 

dt m dr dr dk t 



-;£/^'<W-W>A.£/fer, t) . P> 

Accordingly, we expect that gradient corrections result- 
ing within the standard approaches to the BE are just 
the virial corrections. Our plan is to derive gradient 
corrections from the quasiclassical limit of the quantum- 
statistical transport equation and rearrange them into 
the last term of Eq. (j|) . 

Gradient corrections to scattering integrals have been 
already studied by many authors, the most extensive 
studies were devoted to the effect of electric field on colli- 
sions. Starting from BarkerLa, the field effect on scatter- 
ing has been discussed for number of scattering mecha- 
nisms within various approaches which include the super- 
operatorjarojection technique, Refs. the Lcvinson 

equatiorJla derived from the BBGKY hierarchy, Ref. [l9| , 
and non-equilibrium Green's functions, Refs. |2C|-|25|. 

It is important to note that these different approaches 
resulted in contradictory predictions. Let us focus on 
the (linear) gradient correction to the scattering integral 
which we claim to be equal to the virial correction [the 
last term in the r.h.s of (0)]. We will assume only self- 
energies in the Born or Puff-Whitefield approximation so 
that the virial corrections are absent, i.e., we expect the 
gradient corrections to be absent, too. Surprisingly, this 
is not true about gradient corrections one finds in print. 
One can split all papers into three groups: 

1. In Refs. [16| , |Tt| , p~9|j2^ . linear gradient corrections 
were found. 

2. In Ref. ^5] linear gradient corrections were found 
but with opposite sign than in the first group. 

3. In Refs. |20|-f23| no linear gradient corrections were 
found. 

Unless we will understand why the gradient corrections 
appear for the Born or Puff-Whitefield approximation 
and why there are contradictory predictions from alter- 
native approaches, we can hardly use the gradient cor- 
rections as a start line of our approach to the virial cor- 
rections. 

To identify the origin of the three contradictory predic- 
tions, let us inspect what kind of the transport equation 
is specific for each group. Although various approaches 
have been used by authors, we will describe all the three 
groups within a common dialect of Green's functions. 

1. Integro- differential transport equation for the reduced 
density matrix 

The integro-differential transport equation for the re- 
duced density matrix is obtained— iby the Generalized 
Kadanoff and Baym (GKB) ansatzcJ implemented in the 
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time diagonal of the, integro-differential Kadanoff and 
Baym (KB) equationtl. ^From the GKB ansatz, time ar- 
gument t of the reduced density matrix results retarded 
compared to the time argument T of the scattering in- 
tegral, i.e. t < T. In the integro-differential equation, it 
is natural to identify the instant of the collision with T. 
Using the linear expansion p(t) = p{T) + ^(t — T) one 
then obtains a Boltzmann-like scattering integral from 
p(T) and gradient corrections from jjjr(t — T) . 



2. Integral transport equation for the reduced density matrix 

The integral transport equation for the reduced den- 
sity matrix is obtained by the GKB ansatz implemented 
in the time diagonal of the integral GKB equationcS. In 
the integral equation one has a freedom to identify the 
instant of the collision with t and gradient corrections 
emerge from matching of the scattering integral with the 
subsequent propagation. The opposite signs found in 
groups 1 and 2 thus follow from the fact that in group 1 
authors extrapolate along the initial state of the collision 
while in group 2 authors extrapolate along the final one. 



3. Integro-differential transport equation for the quasiparticle 
distribution 

The transport equation for the quasipaxticle distribu- 
tion is obtained by the original KB ansatza implemented 
to the KB equation at the quasiparticle pole. From the 
KB ansatz, the time argument of the quasiparticle dis- 
tribution results identical to the time argument of the 
scattering integral. This time is naturally identified as 
an instant of the collision event and no gradient correc- 
tions appear. 

The contradictory predictions of the gradient correc- 
tions were found to follow, from different treatments of 
quasiparticle corrections. Ell As shown in Ref. the lin- 
ear gradient correction to the rate of scattering out is 
nothing but the time dependence of the wave-function 
renormalization. In other words, gradient and quasipar- 
ticle corrections are linked together. The approach that 
takes care of quasiparticle features and naturally leads 
to the Landau-Boltzmann type of transport equation for 
quasiparticles is the one of the group 3, i.e., based on 
the KB ansatz. The absence of the gradient corrections 
for the Born approximation is also in agreement with our 
expectation. 



D. Non-virial gradient corrections 

Among studies of the gradient corrections that are not 
devoted to the field effect we want to point out the pa- 
per by Kolomiets and Plyuikcca. In the quasiclassical 
limit, they have evaluated the scattering integral from 



the self-energy in the second-order approximation of the 
electron-electron interaction keeping all gradient terms. 
Similarly to the spirit of classical virial corrections, they 
have expressed gradient corrections in terms of effective 
shifts in space and momentum. The scattering integral 
they have derived thus reminds the one we are looking 
for. 

In spite of formal similarity, the corrections derived 
by Kolomiets and Plyuiko are not the virial corrections. 
Kolomiets and Plyuiko have used a mixed approach im- 
plementing the KB ansatz in the time diagonal of the 
integro-differential form of the GKB equation. These two 
steps are not compatible as the KB ansatz includes the 
quasiparticle distribution while the time diagonal (or in- 
tegral over energy in Wigner's representation) provides 
drift terms for the reduced density matrix. The gradient 
contributions they have found are thus a strange form of 
the quasiparticle corrections. 

Again we came to the same moral. If we want to put 
equality between the virial and the gradient corrections 
to the scattering integral, we have to work fully in the 
quasiparticle picture. 



E. Presented approach 

Our approach is based on two limits applied to non- 
equilibrium Green's functions. The first one.is the quasi- 
classical limit that is inevitable for the BEH.The second 
one is the limit of small scattering rateso This limit 
restricts validity of quasiparticle corrections to weakly 
renormalized systems. The validity of virial corrections 
is restricted to the second virial coefficient. 

The limit of small scattering rates from the quasiclassi- 
cal transport equation for non-equilibrium Green's func- 

In 



tion has been already presented in detail in Ref. 14 



fact, the transport equation for quasiparticles, Eq. (70) of 
Ref. [l4], reduces to the modified BE (Q) . From this point 
of view, this paper is just an application of the method 
derived in Ref. |l4[ On the other hand, in Ref. [l4| au- 
thors discussed only simple self-energies on the level of 
Born approximation which provide quasiparticle correc- 
tions but no virial corrections appear. Accordingly, there 
is no collision delay. Here we use a more complex self- 
energy, averaged T-matrix approximation, that results in 
non-trivial virial corrections similarly as it was done in 
Ref. [l5| for two-particle scattering. 

In spite or because of many similar and closely related 
studies of gradient contributions to transport equations, 
we feel that we have to start our treatment directly from 
non-equilibrium Green's functions instead of recalling al- 
ready achieved results. This is because these treatments 
differ in seemingly tiny details that become important as 
soon as one wants to keep trace of gradients and quasi- 
particle corrections in the same time. 

Now we can specify our aim. The precursor of equation 
(111) is the quasiparticle Boltzmann equation in semicon- 
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(3) 



ductors, Eq. (70) of Ref. [lj], which reads 

df de df de df 
dt dk dr dr dk 

Except for transport vertex <j < , all components are de- 
termined by retarded self-energy a R . First, 7 = —2ima R 
provides inverse lifetime, ^ = zj £ . Second, the quasipar- 
ticle energy is given by e — - — h (f> + Kea R . Note that 
unlike in Paper I, e includes the potential We made 
this change of convention to comply with the convention 
of Ref. If 

It will be easy to show that for homogeneously dis- 
tributed point impurities the velocity simplifies to the 
form in ([!]), J| = z^, and similarly do the force, 
= §f • Apparently, quasiparticle corrections are ex- 
plicitly present in equation (|j|). The viral corrections, if 
present, are hidden in the transport vertex a < . The fo- 
cus of our interest thus will be to show that the transport 
vertex can be rearranged to the form of the scattering-in 
integral from (g), 

- 1 2tt2 



ZG, 



1 2tt_ 2 



dp 
dp 



5(\p\-\k\)f(p,r,t) 



\k\)A t ^f(p,r,t). (4) 



F. Content 

Paper is organized as follows. In Sec. II, we briefly 
review the method of Ref. |l4|. We introduce the non- 
equilibrium Green's functions, the quasiclassical limit, 
and the limit of small scattering rates. With help of these 
tools the Boltzmann-like transport equation is derived. 
In Sec. Ill, the self-energy and the transport vertex are 
specified within the averaged T-matrix approximation. 
The most essential part of our paper is in Sec. IV, where 
we evaluate the scattering integral including its gradient 
corrections. In Sec. V, these gradient corrections are rear- 
ranged into a form that is identical to the one intuitively 
expected and equation (Q) is recovered. In Sec. VI, single- 
electron observables are valuated. Electron density and 
current are discussed in detail. We also discuss the den- 
sity of energy. Sec. VII includes summary. In Appendix 
we discuss relation between virial and quasiparticle cor- 
rections. 



II. TRANSPORT EQUATION 

Our focus of interest is the transport vertex a < . The 
treatment of the transport vertex, however, is intimately 
connected to the transport equation itself. Therefore we 
find it profitable to briefly review derivation of equation 
(||) . This also gives us a room to introduce a necessary set 
of equations for non-equilibrium Green's functions and 
their relations to components of equation (||). 



A. Generalized Kadanoff and Baym equation 

Our starting-point is the generalized Kadanoff and 
Baym equations 



G R E<G A , 



(5) 



that is accompanied by— the Dyson equation for retarded 
(advanced) propagator! 2 ^ 



(Go 



G 



R,A 



(6) 



where the inverse free-particle propagator reads 

<f>(x 1 ,t i ))s(l-2). 



(7) 



Here, numbers are cumulative variables 1 = (ti,x\), 
time and space. Matrix products include integrations 
over time and space, C — AB means C(l,2) = 
/ dt^dx^A{\ : 3)5(3, 2). The e is the free-electron kinetic 
energy and <f> is a potential. Our sign convention for the 
correlation function G < is 



G<(l ) 2) = Tr(^ t (2)^(l)) 



(8) 



where p is the grand-canonical averaging operator, and 
tp^ and tp are creation and anihilation operators, respec- 
tively. 

Equations (|^-||) are general identities and definitions. 
The particular physical content of these equations is spec- 
ified by self-energy T, R ' A and transport vertex S < . The 
transport equation (||) follows from set with no re- 
gards to a particular form of self-energy, except that the 
scattering rate connected to T, R - A and is supposed 
to be small in a sense specified below. Moreover, per- 
turbations in the system of electrons are supposed to be 
smooth in time and space so that corresponding gradi- 
ents are also small. Thus we leave specification of the 
self-energy for next section, and focus on the quasiclas- 
sical limit (small gradients) and limit of small scattering 
rates of equations 



B. Wigner's representation 

In the quasiclassical limit, all operators are conve- 
niently described in Wigner's mixed representation 

a{oj,k,r,t) = J dTdye tUT - iky 

x A(t + 7-/2, r + y/2, t - r/2, r - y/2). (9) 

We use convention that lowercase denotes Wigner's 
transform of operators denoted by uppercase. This con- 
vention is in agreement with the convention used in Pa- 
per I. Indeed, site-diagonal operators used in Paper I (like 
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^R,A,<^ y an j rpR,A^ m Wigner's representation are mo- 
mentum independent and equal to their local elements 
o- R ' A >< , v and t R,A introduced in Paper I. 

The transformation (|J) mixes together left and right 
arguments of the function A, therefore it complicates ma- 
trix products. Keeping only the first gradients in time 
and coordinates, the matrix product C — AB in the 
mixed representation reads 



c = ab + — [a, 6], 



(10) 



where the rectangular brackets denotes Poisson's brack- 
ets 



. ,. da db da db da db da db 

' ' ' ' duj dt dt duj dk dr dr dk 



(11) 



C. Propagation of single quasiparticle 

In the quasiclassical limit, one restricts the assumed 
fields <f> to those that vary slowly in time and space. In 
this case, the field is expected to have a pronounced effect 
only on long trajectories of electrons. Such long trajec- 
tories necessarily include number of collisions with impu- 
rities. A propagation between subsequent collisions and 
individual collisions themselves happen on a very small 
time and space scale, thus they should be nearly the same 
as in a homogeneous and stationary potential, i.e., in the 
absence of the field. There might be also a small effect of 
the field on these microscopic scales, however, this field 
effect can be handled as a correction linear in gradients 
of the field. 

Propagators G R ' A describe a motion between colli- 
sions and also internal dynamics of collisions. Thus 
we have to find them up to linear order of gradients. 
From the Dyson equation (^) and its alternative form 
G r a (Gq 1 — S^rn 4 ) = 1 one finds that propagators are 
free of gradientsa 



<J 



R,A 



1 



UJ — £ — (j) — G 

Being a complex function, 



R.A ■ 



a 



R.A 



(12) 



(13) 



the self-energy describes two kinds of phenomena. Its 
imaginary part ^7 descries scattering out of the state of 
momentum k. Its real part a renormalizes energy of the 
single-particle-like state. The renormalized (quasiparti- 
cle) energy is given by a position of the pole of G RA at 
the real axis 



■<j> + a{e). 



(14) 



This quasiparticle energy is an ingredient of the quasi- 
particle BE (|). 

In the limit of small scattering rates (small 7 expan- 



sion, 



7 



0), the pole of the propagator sits 



close to the real axis. Then one approximates the spec- 
tral function 



-2Inig J 



by its limiting value 

a = 2ttzS(uj — e) + 7 Re 



1 



(u - e - (f> - iO) 2 



(15) 



(16) 



Within the quasiparticle picture, in the spectral func- 
tion a(u>, k, r, t), the 5 function represents a singular con- 
tribution from quasiparticle state of momentum k (at 
point r and time t). The second term describes vari- 
ous projections of other quasiparticle states with energy 
e(jp, r, t) — lu into momentum k, in other words, the sec- 
ond term is an off-pole contribution. Within the quasi- 
particle picture, these two parts of the spectral function 
has to be treated separately. 

The norm z of the <5 function is the wave-function 
renormalization 



1 



dcr(uj) 



(17) 



D. Completed collisions 

The generalized Kadanoff and Baym equation (||) in- 
cludes all phenomena we are interested in, but in a hidden 
form. The simple link between Green's function formulae 
and intuitive approach to the BE with virial corrections 
is obscured by the fact that these two approaches deal 
with different objects. While the correlation function 
G < includes both, the pole and off-pole, contributions, 
the quasiparticle distribution is related only to the pole 
part. To recover the BE equation, one has to separate 
these two parts. 

Within Green's function the pole and off-pole propa- 
gation are described by a single transport equation (|^). 
In contrast, a diffusion of particles described by the LHS 
of the BE applies only to the pole part, while the off- 
pole part seems to be missing. In fact, it is not missing. 
The off-pole propagation is hidden in scattering integrals 
and relation between observables and the quasiparticle 
distribution. 

The possibility to move the off-pole propagation into 
the scattering integral follows from a hierarchy of time 
and space scales in the system that are inevitable for the 
theory of Boltzmann type. Diffusion of quasiparticles is 
well defined on a hydrodynamical scale that includes a 
large number of impurities and a large number of colli- 
sions per particle. On the hydrodynamical scale, there 
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is an appreciable effect of the field <j> which may cause a 
transfer of quasiparticles over long distances. Such a mas- 
sive changes of the system state is effectively described 
with help of a differential transport equation that bal- 
ances drift with dissipation. In contrast, individual colli- 
sions happen on a microscopic (local) scale on which the 
effect of the field <\> is small so that it can either be ne- 
glected or included by corrections linear in its gradients. 
Thus, a subdynamics on the microscopic scale can be 
integrated through and approximated by effective scat- 
tering rates. 

Clearly, to recover the BE we have to separate hy- 
drodynamical and microscopic scales. This separation is 
equivalent to the separation of the pole and off-pole parts 
based on a non-equilibrium,papdification of the expansion 
in small scattering rates. EalJ Alternatively, the separa- 
tion of pole and off-pole parts of G < can be based on the 
idea of completed collisions. In this paper we focus on 
the microscopic mechanism of collisions, thus we follow 
the completed collision approach. 

For the purpose of motivation, we assume for a while 
a homogeneous system that is diagonal in momentum 
representation. Let us take a look on a time-diagonal el- 
ement of G < , say G < (t,t;k). In the transport equation 
(|l7|), the transport vertex represents the last collision 
due to which a wave function of electron gained a com- 
ponent of the momentum k. This component can belong 
either to an asymptotic state that will form new effective 
quasiparticle state or to the off-pole contribution due to 
some other state p. The asymptotic state is on the energy 
shell, thus it will live on the time scale of the quasipar- 
ticle life time. The off-pole contribution will vanish on 
the scale of a quasiparticle formation time. The latter is 
much shorter than the former. To distinguish whether 
the contribution to G < (t, t;k) is of pole or off-pole na- 
ture, one can monitor a vicinity of the time t, including 
a close future, to figure out which part will survive and 
which part soon disappears. This procedure corresponds 
the approach of Fermi golden rule, where one also inte- 
grates through a collision into future and uses matching 
of asymptotic states to identify a product of the com- 
pleted collision. 

In accordance with the causality principle, all time in- 
tegrals in ([It]) run only over the left part of the time 
axis, i.e., for times smaller than t. To monitor a close 
time vicinity, we rearrange the transport equation ( p"^ ) 
as 

G< = i(G fl - G A )Y> < G A ~ ^G R Z<(G R - G A ) 
+ ^G R Z < G R + ^G A I1 < G A . (18) 



In the added term 

S< = -G R ^ < G R + -G A ^ < G A , 
2 2 

the time integration runs into future, 



(19) 



2~ < (t ll t 2 )= I dt' I diG R (t 1 ,t>)Z<(t',t)G H (i,t 2 ) 

poo 

dt / dt'G A {ti,t')Y l < {t' ',i)G A (i,t 2 ). 
> Jti 

(20) 

The same integration into future also appears in the other 
contribution to G < 



A < 



1 



(G 1 



G /1 )S < G 



<r* A 



1 



G £ (G 



(21) 



On the time diagonal, — t, the time integrations 
over t' and t in S < do not overlap, see (20). The time 
scale of the integration is determined by the time scale 
of S < that can be identified with the quasiparticle for- 
mation time.E3 Thus S < is dominated by the short-time 
(off-pole) contributions. 

The time integration in A < extends into the future in 
a way that reminds Fourier transformation to energy of 
the asymptotic state, 



G R (t,t';k)-G A (t,t';k) 



-ic 



-ie k (t-t') 



(22) 



Thus A < is dominated by the long-time (pole) contribu- 
tions. 



E. Pole and off-pole contributions 

Splitting the correlation function G < into A< and S < 
is an ideal starting point to a non-equilibrium modifi- 
cation of the expansion in small scattering rates. This 
can be seen in equilibrium, where one can use the energy 
representation 

1 



\ < (u,k) = f FD (Lu)-a(ui,k) 

-> f(k)z(k)2nS(u - e). 



(23) 



We have used <J < (uj) = fpD{^)j(uj), and The arrow 
shows a value in the limit of small scattering rates. 

Note that X K = fFD^a 2 ^, while g < = fFD a - In the 
limit of small scattering rate, 7 — > 0, the spectral function 
a approaches 5 function as Lorentzian, while s = \a 2 ^ 
approaches i5 function faster, 



a(w, k) = 
s(u>, k) = 



7 



2J_ 

((w - e - (7)2 + I 7 a) 



2 ' 



(24) 
(25) 



In the off-pole region \w — e| 3> 7, the spectral function 
a has a tail linear in 7. In the limit of small scattering 
rates, 7 — > 0, this tail results in the off-pole correction, 
the second term in (|l6|). The function s in the off-pole 
region is proportional to 7 3 , thus its limit is a pure 5 
function without the off-pole term. 



G 



The function contains the off-pole part. In equilib- 



rium, 



£<(w,fc) = / F . D (w)7(w)Re(«-e-ff fl ) S 
-> /fd (w)7(w)Re (w - e + iO) 2 . 



(26) 



Due to the off-pole nature of this contribution, the Fermi- 
Dirac distribution in (26) cannot be associated with oc- 
cupation of the state k. 

The comparison with equilibrium shows that the quasi- 
particle distribution relates to A < while has to be con- 
structed indirectly. In the spirit of the BE, we will treat 
A < within a differential transport equation while ^ < will 
be turned into a local functional. 

Equilibrium relations ( ^3^2^ ) can be easily generalized 
to quasiclassical limit. As we have shown, the spectral 
function a is free of gradients, see (|l3| ) and (|l5|). The 
function S = A + \G R TG R + \G A TG A is also free of 
gradients since gradient expansion of symmetric terms 
like G R TG R has no gradients. Therefore, similarly to 
the spectral function a, the function s just follows energy 
bottom defined by the field 4> 



1^3 
2 ' 



((c-e-^-rja + i^)' 



(27) 



In the limit of small scattering rate, 7 — > 0, the function 
s reduces to the first term of the spectral function, 



s = 2t: zS(lo — e). 



(28) 



In equilibrium, the pole part A < is proportional to the 
function s. Out of equilibrium, we can expect similar 
behaviour and introduce local distribution as 

A< (w, k, r, t) = f%(uj, k, r, t)s(w, k, r, t). (29) 

In the limit of small scattering rates, the function s turns 
to the 5 function and one can eliminate the energy argu- 
ment of the local distribution. In this way we can define 
the quasiparticle distribution as the pole of the local dis- 
tribution. Briefly, in the limit of small scattering rates 
the pole part of the correlation function reads 



A < (w, k, r, t) = f(k, r, t)2nzS(uj — e). 



(30) 



The non-equilibrium generalization of the off-pole part 



F. Boltzmann equation 

The BE is recovered from (|2~l|). First, we turn ( pl| ) into 
differential form multiplying it by G^ 1 from the l.h.s. and 
by G^ 1 from the r.h.s. and subtracting the two forms 



z(G^A< 



A<G^) 



= - (E<A + AT.< - TG A T,<G A - G R T,<G R T). (33) 

By the gradient expansion, this equation simplifies as 

[u - e - <j> - a, A<] - ~[ 7 , aga<] = a<s - 7 A<, (34) 

where g — Keg R . In the limit of small scattering rates, 
the second term on the l.h.s. vanishes. After a substitu- 
tion of (|30"|), equation (|i]) turns into equation (||). For 
details see Refs. |U||. 

To deal with the transport equation (^), we need: life 
time t, quasiparticle energy e, wave-function rcnormal- 
ization z, and transport vertex tr < . Since all these func- 
tions are related to the self-energy, further progress re- 
quires to specify the self-energy. For elastic scattering on 
impurities, the retarded (advanced) self-energy depends 
only on the retarded (advanced) propagator G R . The 
transport vertex depends on both propagators, G R and 
G A , and on the correlation function G < . Propagators are 
given by (|l3|). The correlation function has to be decom- 
posed into two parts according (|l^ ) which in Wigner's 
representation reads 
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A< + r 



(35) 



The pole part A < is linked to the quasiparticle distri- 
bution via (|30|), the off-pole part is self-consistently 
evaluated from the transport vertex a < via (p2|). 



III. AVERAGED T-MATRIX APPROXIMATION 

Here we specify the self-energy. As in Paper I, we as- 
sume non-interacting electrons scattered by neutral im- 
purities. The impurity acts on electrons by a jsingle-site 
potential introduced by Koster and SlaterEilO. Individ- 
ual scattering events are treated within the T-matrix. 
Formulas for the T-matrix in homogeneous systems are 



follows directly from its definition (ng). Since sym- quite common,EI our focus will be on gradient contribu- 



metric terms has no gradient contributions 

i < =o < \{g\^g\). 

In the limit of small scattering rates 

1 



= er Re 



{uj-e-(j)-a + iO) 2 



(31) 



(32) 



tions in inhomogeneous systems. 



A. Retarded self-energy 

In the Koster-Slater model, an impurity at position r is 
characterized by a potential restricted to a single orbital 
\r) at site r, 



V r = \r)v{r\. 



(36) 
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Corresponding retarded T-matrix reads 
T R = V r + V r G R T R . 



(37) 



Iterating (|37|) one can see that the T-matrix is also single- 
site function, 



T 1 



V) 



1 - (r\G R \r) 



(r\ 



(38) 



The T-matrix does not depend on difference coordi- 
nate, therefore its mixed representation relates only to 
double-time structure 

|r>t*(w,t,r)<r|= / dre i ^\r)t R (t + ~ t - I, r)(r|. 



(39) 

Similarly, the local element of the propagator also de- 
pends only on energy, not on momentum. In Wigner's 
representation, the local element of the retarded prop- 
agator {r\G R {t u t 2 )\r) = G R {r,tx,r,t 2 ) = G R (r,h,t 2 ) 
transforms into 



dk 



■g R {uj,k,r,t). 



The T-matrix in mixed representation reads 

R . . v 

t (w,r,t) = tt • 

1 - vg H {uj,r,t) 



(40) 



(41) 



There are no gradient contributions. This can be checked 
directly by explicit gradient expansion of (|38|). 

The retarded self-energy is defined as a mean value of 
the T-matrix 



S fl = J drc(r)T r R , 



(42) 



where c(r) is a concentration of impurity per site on 
site r. This approximation is called the self-consistent 
averaged T-matrix approximation (ATA). Unlike in Pa- 
per I, we do not use the subscript self here since the 
self-consistent form is a natural starting point in trans- 
port theory. Non-sclf-consistent values are introduced 
and denoted below. 

In Wigner's representation (E|) reads 



a R {tu,r,t) =c{r)t R (LO,r,t). 



(43) 



which in Wigner's representation reads 
o< = ct R g<t A 
+ c l -g<[t R ,t A ] 

+ c l - S<] 



Here. 



dk 
(2^V 



■G K {uj,k,r, t) 



(45) 



(46) 



is the local element of the correlation function. The Pois- 
son bracket used in ( f45| ) as a short-hand for gradient cor- 
rections in general includes space derivatives combined 
with derivatives with respect to momentum, see ( pT| ) . In 
([45|), however, none of functions depends on momentum 
so that the only gradient contributions come from time- 
derivatives. 

The transport vertex a < has three basic components 
that have distinguishable physical content. First, there 
is a non-gradient term 



a< = ct R g<t A . 



Second, there is the term 



a<=ci~g<[t R ,t A ] 



(47) 



(48) 



which formally brings gradient corrections to the scatter- 
ing rates. Below we show that this term vanishes. Third, 
there are two complex conjugate terms 



(49) 



which contribute if the quasiparticle distribution / has 
non-zero gradients. This last term results in non-local 
corrections of the BE (|l|). 

Now the set of equations is complete. 



IV. SCATTERING INTEGRALS 

The scattering-out integral z^/ £ f and the non-gradient 
part of the scattering-in integral zaf are dominant. They 
are the only non-zero terms of the BE in the absence of 
the perturbing field <j) and their balance determines how 
the BE is capable to describe equilibrium. 



B. Transport vertex 

The transport vertex a < in the self-consistent ATA 
depends on the correlation function G < as 

£< = / drc{r)T R G<T A , (44) 



A. Non-gradient part of the scattering integral 

In the quasiclassical limit, the frequency of the per- 
turbing field is much smaller than the relaxation time of 
the system, i.e., all perturbations are on the long-time 
scale. Therefore, the consistency of z^/ £ f and zaf is cru- 
cial and it has to be checked in very detail. 
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On the level of non-equilibrium Green's functions, the 
scattering in and out are consistent if the self-energy 
is given by the T-matrix that obeys the optical theo- 
rem. The ATA obeys it. The consistency on the level 
of Green's functions does not imply a consistency within 
the BE. The scattcring-in integral in the BE include addi- 
tional approximations (|3(]) and (j32f), while the scattering 
out is evaluated without them. Since (|3(]) and (|32] ) follow 
from the limit of small scattering rates, we have to make 
corresponding approximation for the scattering out. 



1. Scattering out 

The limit of small scattering rates, 7 — > 0, is conve- 
niently discussed in terms of the non-self-consistent ATA. 
Indeed, sending a R — > in propagators g R of the self- 
consistent T-matrix t R turns into the non-self-consistent 
one t§ Q . 

^From ( p2[ ) one can see that 

„R(, ,\ 1, , a ~R\ 



'M=0oo(«-*-O> 



(50) 



where subscript 00 denotes no field and no self-energy in 
propagator, [in Paper I, these non-sclf-consistcnt func- 
tions are without subscript] 



0oo M = 



1 



uj — e + iO 



(51) 



The self-energy a R is then expressed in terms of the non- 
self-consistent self-energy 



a R - rf R - r 
°oo — C( oo — c 



1 - Wo 



(52) 



as 



da 
du> 



da, 



on 



dui 

which is equivalent to 
1 



uj—e—<j)~a 



1 - 



duj 



dciQo 



1 



duj I uj—e 



du 



(57) 



(58) 



Accordingly, we have found that the self-consistent and 
non-self-consistent approximations are linked via the 
wave-function rcnormalization 



- = zj s = 700(e) = c(-2)Imt$,(e). 



(59) 



The non-self-consistent T-matrix satisfies the optical 



theorem lmt R 



t R (lmg R Q )tQ Q , which easily follows from 



([y]) and the complex conjugacy of retarded and advanced 
functions, = {t R )* . The scattering-out rate thus can 
be expressed in form of a sum over individual scattering 
rates into all accessible finite states 



dp 



■2nS(e k - e p ). 



(60) 



The 5 function in the r.h.s. results from —2^300(0;) = 
2ir8(w — e). From (|60|) one identifies the scattering rates 



Ppk = c \t R {e k )\ 27r5(e fe -e p ). 



(61) 



Since there are no gradient contributions to the scatter- 
ing out, the quantum-mechanical scattering-out rate ( |6l| ) 
is of the same form as the intuitively expected classical 
scattering-out integral in the modified BE fah. 



2. Scattering in 



a R {L0)=a R {u-4>-a+- 1 ). (53) 



In the limit of small scattering rates we linearize in 7, 

■ (54) 



a R {u>) = „ R i"-<t>-v) + h da ™ 



uj — m—a 



2 ' duj 

The imaginary part of (p4), 7 = — 21ma R , reads 



7oo(w -4>-<j). (55) 



U) — <j) — (J 



At the pole value aj~e~e + cj) + <r(e), the argument 
of 700 simplifies 



(1 • ^) _ = 7oo(e). 



(56) 



The factor 1 + ^f 22 - 1 is iust the wave-function renor- 
malization z, see Appendix A of Paper I. Indeed, from 
real part of one finds 



Now we evaluate the non-gradient part a< of the 



scattering- in integral (|47]). We show that it results ex- 
actly in the scattering-in inte gra l from the BE (|l|) with 
the scattering rates given by (pl[). In this way, the con- 
sistency of scattering in and out will be checked. 

Substituting g K from ([l6]) with g < = X < + from 
<M) and @, one gets 



c\t R \ 2 r 



= c\f 



+ c r 



dp 
dp 



zf(p)2nS{uj - e p ) 



(2tt) 3 ng (uj-e p -4>-<T + i0y 



(62) 



The first term follows from A < and is the dominant one. 
The second one is the off-pole correction due to . 

The T-matrix in (|62|) includes propagators with poles 
out of the real energy axis, shifted by 7/2. In the limit 
of small scattering rates, it is advantageous to use the 
non-self-consistent ATA as the reference point 
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i fl (w)=tg,(w-0- 



77) 



du> 



(63) 



^From this approximation one finds that the square of 
the absolute value at the pole reads 



\t R (e)\ 



1 r)t A 
^00 duj 



(64) 



The second term in the bracket in ( |64| ) can be expressed 
in terms of the collision delayL-3 (1-25) 



A, 



1 r)t R 
Im 1 dt °° 



as 



*00 9cj 



|t*(<0| 



1 



(65) 



(66) 



This relation between the self-consistent and non-self- 
consistent T-matrices has a form of virial corrections. 
These are not, however, the virial corrections we are 
looking for. The term At w iH be cancelled by the wave- 
function renormalization and the off-pole contribution. 

Let us denote by <Jq the non-self-consistent-like scat- 
tering vertex 



°oo( £ k) = c\t$b(£ k - <f> -a) I 



dp 



■zf(p)2ir5(e k - e p ) 



c|*oo(efc) 
dp 



dp 



(27T) 



(2tt) 3 

8 *W(p)- 



f(p)2nS(e k 



(67) 



Apparently, er^ is compatible with the scattering out. 
To prove the compatibility of a^g with f'y, we have to 
show that all corrections (linear in 7) following from the 
non-self-consistency mutually cancel, i.e., za< g = (Jq . 

There are three contributions to the correction za^ g — 
(7 qq : (i) the second term in (|62| ) , (ii) the At correction to 
the square of the T-matrix, and (iii) the wave-function 
renormalization in front of af . Remains to show that the 
sum of these three corrections vanishes. We will neglect 
higher order terms resulting from products of individual 
corrections. 

First, we rearrange the second term in (^). To this 
end we use that o < does not depend on the momentum 
p and move it out of the integral. The integrand can 
be rearranged is spirit of Ward's identities with help of 
energy derivative as 



Re 



1 



{u>-e-(f>-a + iO) 2 



div 1 



da_ 



(68) 



By the integration over p : relation ( pq ) turns into the 
local propagator needed in the second term of (p2[) 



dp 



Re- 



1 



(2tt)3 (e k 



(69) 



The second term of (f32|) thus can be rearranged as 



dp 
(2^ 



cr5. Re 



1 



-za< g c\t^ {e k )\" ( 1 - -^7 



{Ek-e P -<t>- cr(e k ) + iO) 2 
r z 



(70) 



Now we can collect all terms which contribute to the 
non-gradient part of the scattering in 



za nA £ k) = Z 



— \ < 1 \ \±R |2 ^500 



A, 

T 

duj 



za< g (e k ) 



c\t 



R |2 
00 1 



dffoo 
duj 



o-oo( £ fc) 



= o-^o(e fc ). 



(71) 



In the last but one step, we have neglected the cross- 
correction At x c^q!)! 2 -! 22 - and terms quadratic in 



c\t, 



R 12 dgoo 
00 1 du 



which are higher order in the limit of small 
scattering rates. In the last step, we have used the de- 
rived optical theorem (I-B3), [proved also in Appendix of 
this paper, see ( A17[ )] 



At 
r 



dw 



c\t 



R |2 
00 1 



dgoo 
duj 



(72) 



Briefly, we have shown that in the non-gradient scat- 
tering in and out, the wave- function renormalization z 
compensates in a consistent manner. This compensa- 
tion shows that the off-pole part S < of the correlation 
function is capable to rebound consistently the off-pole 
portion of the particle propagation. 

The non-gradient parts of scattering integrals 



zj(s k ,r, t) 
za< (s k ,r, t) 



dk 



(2tt) 3 
dk 

W) 



P kp f{k,r,t), 
3 P P kf(p,r,t), 



(73) 
(74) 



have scattering rates identical to their non-self-consistent 
counterparts 



P P k = c|^ (efc)| 2 27T(5(efc - e p ) 

= \ 2 -^m-\k\). 



(75) 



These non-gradient parts are identical to those used in 
the intuitive BE (0). 
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B. Gradient corrections to the scattering rate 

The gradient term of the self-energy crp given by J4S| ) 
does not include gradients of the quasiparticle distribu- 
tion. Accordingly, dp can bring corrections to scattering 
rates P p k only. In reality o-p = 0, i.e., no correction to 
scattering rates appears. 

For impurity potentials independent from time, v = 
const, the Poisson bracket in (J48) equals zero. Indeed, 
the both T-matrices depend on the time t only via the 
internal potential <fi which always enters the T-matrix in 
the form of u> — <fi, i.e., t R ' (lj, t) — t R > (lo — <fi). Any two 
functions of this property has zero Poisson bracket, 

da db da db 

8uj dt dt dui 
da db d(f) da d<fi db 

duj doj dt duo dt doj 
= 0. (76) 



[a(w — 4>), b(uj 



Accordingly, 



0. 



(77) 



Briefly, amplitudes of scattering rates are not modi- 
fied by gradient corrections for time-independent impu- 
rity potentials. This is in agreement with the intuitive 
expectation used in Paper I. 



The function zS(u> 
[zS(uj — e) = S(ui - 



- e) depends on time only via to — <p 
e — <f> — cr)], therefore according to 



~3T~§Z ) 27Tz<5(u; — e p ) = 0. The gradient 



correction thus depends exclusively on gradients of 
the distribution function, 



A 



= -C\t 



R\2 



dp 



1 dt R d 
2vz5{u J -E p )lTn- [E — -f{p). 



ID 



Finally, we simplify this gradient correction with help of 
(|66|), scattering rate (|6l]) and the collision delay (|65|) as 



dp 
(2^)3 



p, 



pk 



'00 

df 
dt 



^ 1 27rzS(cj-s)A t ^ i f(p) 



(27T) 



A, 



(82) 



In the second line, we have neglected higher order terms 
using the approximation z — A t /t w 1 . 

This gradient correction to the scattering integrals has 
the form of virial corrections, the last term of (||). In 
particular, no other gradient than the time derivative 
of the quasiparticle distribution appears, and this time 
derivative is weighted with the collision delay and the 
scattering rate. 



C. Gradient corrections to scattering in 

The gradient correction er^ brings true non-local con- 
tributions to the scattering-in integral. Via gradients of 
the local correlation function the correction de- 
pends on gradients of the quasiparticle distribution. Here 
we show that from <t^ one recovers virial corrections in 
form discussed in Paper I. 

First, we write (}49|) in explicit form 



R+A 



ct"t 



C [t«[ a Im-^[t«,ff<] 



= -c\t R \ 2 lm^ | ^— — 



t R[ 

1 fdt R d dt R d 



t R V duj dt dt duj< 



9 < - 



(78) 



Now, we substitute the dominant part of the local cor- 
relation function (the off-pole part leads to higher order 
contribution in small 7) into ( |4S| ) 



9 < {u,r,t) 



dp 
(2^ 



■fiP, r,t)2irz5{uj - e) 



(79) 



and interchange derivatives with the momentum integral 



ov = —c\t 



R\2 



dp l_ 

(2tt) 3 t R 



I ' dt R d dt R d 



(80) 



V. RECOVERING THE INTUITIVE 
TRANSPORT EQUATION 

Now we can put together elements of the transport 
equation @ and reconstruct (|l|). 



A. Drift part of the transport equation 

The velocity results from momentum derivative of (|l4| ) 



as 



de de da(u>) 
dk dk duj 
1 



uj—e 

de 



de 
dk 



1 - 



dcr(uj) 



dk' 



(83) 



,2 

For parabolic kinetic energy, e = the quasiparticle 
velocity gains the form used in the BE (Q), 



de 
dk 



(84) 



The force acting on quasiparticle is also found from the 
quasiparticle energy, 



de dd> da da de 
— = — H 1 . 

dr dr dr dui dr 



(85) 
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For homogeneous distribution of impurities, the self- 
energy depends on coordinate r only via potential <f>. 
Since the potential cj) can be viewed as a local shift of 
initial of energies, one finds that the self-energy a relates 
to the self-energy o$ = q in the absence of the field cj) as 



a(uj) = (T0 = o(w - 



(86) 



^From (|8(j) one finds that the second term of the force 
( |S5| ) reads 



da 
dr 



da dcj) 
duo dr 



(87) 



Using ( p7|) in (jS5|), one finds that the force has no renor- 
malization 



de 
dr 



dr ' 



which is the form of the force in the BE (Q). Thus drift 
terms of (||) reduce to drift terms of (|l]). 



B. Scattering integral with virial corrections 

To obtain the expected form ([!]), we use that within 
linear approximation 



a(x) 



dx 



A = a(x + A). 



(89) 



Thus the non-gradient and the gradient scattering-in con- 
tributions can be collected into a compact expression 

zaf = za< g (e k ) + za^(e k ) 



dp 
dp 



P P k[f 



df 



(2tt) 3 
t k 2 J (2tt) 3 



P pk f(p,r,t-A t ) 
dp 



S(\ P \-\k\)f(p,r,t- At). (90) 



One can see that the scattering in has the form expected 
from classical assumptions. 

Substitutin g ([59l) for the scattering out, ( po|) for the 
scattering in, (|84[) for the quasiparticle velocity, and ( pSj ) 
for the accelerating force into the asymptotic equation 
(||), the intuitive modification of the BE (|l|) is recovered 
from quantum statistics. 



VI. OBSERVABLES 

To recover the transport equation (|l|) was our major 
task. With respect to applications, one has also to find 
relation between quasiparticle distribution / and observ- 
ables. As already shown in Paper I, these relations in- 
clude quasiparticle and virial corrections. In Paper I we 



have discussed only those observables that can be iden- 
tified from the transport equation via conservation laws. 
Here we extend our treatment to a general single-particle 
observable. 

All single-electron observables can be expressed in 
terms of the reduced density matrix (Wigner's distribu- 
tion function) 



p(k,r, t) 



2^ 



G<(cj,k,r,t). 



(91) 



^From the decomposition 
are given by (|30| ) and ([ 
reduced density matrix as 



/ < = A < + where A < and 
12J), respectively, one finds the 



P = zf - 



dw 



P 



da < 



2n uj — e du> 



(92) 



Note that the transport vertex a K in (^2|) does not en- 
ters the reduced density only by its pole value af but the 
full energy dependence has to be maintained. Since the 
second term is already an off-pole correction, the correla- 
tion function of the self-energy is used only in its lowest 
approximation a < = a^ Q . 

The formula ( |92|) has no explicit gradient terms, how- 
ever, there are gradient contributions hidden in the trans- 
port vertex ct < . There is a question whether one should 
keep these gradient corrections in the off-pole part of for- 
mula (|9|) or not. A general answer is not clear to us. 
On the other hand, with respect to conservation laws 
that test a consistency of observables with the transport 
equation ([!]) , the gradient contributions can be neglected. 
The transport equation does not provide observables but 
only their time or space derivatives, see e.g. the equa- 
tion of continuity (I-C2), therefore any gradient contribu- 
tion to observables enter the conservation law via second 
derivatives that are neglected within the quasiclassical 
limit. Accordingly, we neglect gradient corrections to a < 
in (M. 



A. Local density of electrons 

In Paper I we have derived the local density of elec- 
trons n from the transport equation, see (1-5). For scat- 
tering on the Koster-Slater impurities it was found that 
the correlated density n corr is determined by the ratio of 
the collision delay to the lifetime, see (1-29). Here we 
recover (1-29) directly from (|92"|). 

The local density of quasiparticles reads 



dk 



/(*). 



(93) 



The local density of electrons is given by the integral from 
p over momentum 



dk 



(94) 
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/.From ([ 



and ( |92| ) one finds 
dk 



rf(k) 



(2tt) 3 
dk daoo 



(2tt) 3 8uj 
dk f diu 



(2tt) 3 



f(k) 



27T W — Eft cta-I 



(95) 



where we have used (f58|) for the wave-function renormal- 
ization z. The first term is the quasiparticle density nf roo , 
the second and the third terms are the correlated density 



dk dcr Q 



(2tt) 3 du 
dk f duj 



f(k) 
P da< (uj) 



(96) 



(27r) 3 J 2tt uj — Sk doj 
In the second term we perform the integration by parts 
dp „ 1 



= -Re- 
cta; lu — e (a; — e + zO) 2 ' 

and substitute for the self-energy a^ from (t 



(97) 



dk 



f(k) 



da o 



(2tt) 3 
dk dco 



du> 



-Re 



UJ=Sk 

1 



(2tt) 3 2vr (u - e k + iOJ 



- c Vo^ 



dp 



■zf(p)2n5{uj - £ p ). 



(98) 



Due to the energy-conserving 5 function, the energy u> 
can be easily integrated out. The wave- function renor- 
malization z under the p integral can be omitted as higher 
order in the limit of small scattering rates, 



dk da o 



(27T) 
dk 



dp 



(2tt) 3 (2tt) 3 



f(p)Re 



(e P - £fc + «0) 2 



koo( e p) I 



(99) 



We have used (|1J) to simplify energy arguments. Now we 
can integrate over momentum k using non-self-consistent 
version of 



dk da Q0 
— f(k) 



(27T) 

dp 



-J(P) 



du 
dgoa 



8lo 



|^oo( e p) | 



(100) 



The two terms can be joint. In the first term we rename 
the integration variable k to p, so that both terms will 
have the same name of the momentum argument of the 
distribution function /, 



dp , / da o 
— f(p) 1 



(27T) 



\ duj 



c \ t. 



1 2 dg { 



(10 



00 



dui 



(101) 



Finally, we apply the derived optical theorem ( |72] ) to re- 
cover the relation (1-29) 



dp , > A t 



(27T) 



(102) 



The total density n = rtf rf o + rt corr resulting from the 
reduced density matrix (|92| ) is thus identical to the one 
obtained from the transport equation via the equation of 
continuity, 



dp 



A, 



(103) 



which is Eq (1-48). Briefly, with respect to the electron 
density, the approximative functional for the reduced 
density matrix (92) is consistent with approximations in 
the transport equation (0). 



B. Local density of current 

The particle current is one of the quantities most often 
evaluated from the BE. Here we show that for the Koster- 
Slater impurities there are no explicit virial corrections 
in the functional for current. 

A general formula for the current is 



3 



dk de 
-^3 W : 



dk k 



(27r) 3 dfc r J (27t) 3 to' 
Now we substitute for p from (|92| ) 
dk k 



(104) 



{2n) 3 m Zf 
dk du> k 



-Re- 



a<(u) 



(2tt) 3 2tt to {uj - e k + iO) 2 

dk k / , , 

\2tt) 6 to 



The second term in the first line is zero because its inte- 
grand is an odd function of momentum k. The current is 
thus identical to (1-41) found from the equation of conti- 
nuity. 



C. Local density of energy 

Energy of the system is not a single-electron observ- 
able. Although electrons do not mutually interact, the 
energy of the system cannot be evaluated from the re- 
duced density matrix p. This is because electrons are 
correlated with impurities of unknown positions. Simi- 
larly to interacting systems, the energy has to be evalu- 
ated directly from the correlation function g < , 
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E = 



dto dk 
2^ (2tt) : 



■ojg < (u, k). 



(106) 



To express the energy in terms of quasiparticle distribu- 
tion, we use g < = A< + with X < from (|3^) and ^ < 
from (|3S 



, duj dk ( ( da \ 1 „ , , . 



f(p)2ir6(uj - e p ) 



x Re 



1 



(w - e k + iO) 1 



(107) 



Now we integrate out the energy u) and interchange the 
order of integrations over momentum in the off-shell 
term, 



E = 



+ c 



dk 
(2tt)3 
dp 



e k f(k) 



dk da 00 

— £kf{k) a7 



(27T) 



(2tt) 3 



r Re 



(2tt) 3 (e p - e k + iOf 



(108) 



Similarly as in the case of the electron density, we join 
the second and the third terms, evaluate the integration 
over momentum k, and apply the derived optical theorem 
((72|) to express the local energy in terms of the collision 
delay 



E 



dk , J A t 



(27T) 



(109) 



[In notation of Paper I, Eq. (1-42), the quasiparticle en- 
ergy does not include the potential <j> that is in (1-42) 
explicitly added.] Thermodynamical consistency of the 
energy conservation and correlated density is shown in 
Appendix C of Paper I. 

The BE (JT|) with subsidiary relations fli03| ), ( $Mj) and 
( |109| ) form a basic set of equations that cover most of 
traditional applications of the BE. In all these equations, 
the virial corrections can be included with help of colli- 
sion delay. We remind that this simplicity follows in part 
from the simplicity of the scattering by Koster-Slater im- 
purities. 



VII. SUMMARY 

The intuitive modification of the BE, Eq. ([[]), has 
been recovered from non-equilibrium statistics. To this 
end we have employed non-equilibrium Green's functions 
within which we made the quasiclassical limit and the 
limit of small scattering rates. These two limits are fully 
sufficient, i.e., no unjustified approximations need to be 
made. The non-local form of the scattering integral in 



the intuitive BE has been obtained by unification of non- 
gradient and gradient contributions. 

Single-electron observables as functionals of the quasi- 
particle distribution are provided by the reduced density 
matrix that in the limit of small scattering rates has form 
( p2] ) . It was shown that (^2|) is consistent with the trans- 
port equation (|l|) leading to the correct equation of conti- 
nuity discussed already in Paper I. The density of energy, 
which does not belong to single-electron observables, has 
been treated separately. 

Presented theory has four general features that can 
be transferred to more general models. First, one needs 
sufficiently complex self-energy, the recommended one is 
based on the T-matrix which guarantees number of iden- 
tities related to the optical theorem. Second, for small 
scattering rate, one can use the procedure of Refs. [Ti|j29| 
to derive a quasiclassical transport equation for quasipar- 
ticles. Resulting transport equation includes the quasi- 
particle and the virial corrections. The virial corrections 
are, however, in a form of gradient contributions to the 
scattering integral. Third, the virial corrections are re- 
arranged to the semi-classical form when one recollects 
non-gradient and gradient terms on the scattering inte- 
gral using logarithmic derivatives. Fourth, all logarith- 
mic derivatives should be defined from the T-matrix, i.e., 
from the scattering phase shift. These logarithmic deriva- 
tives have natural interpretations like the collision delay 
discussed here. Authors of this paper have implemented 
presented approach to nuclear matter, see Ref. |35[ 

On the other hand, the discussed scattering on neu- 
tral impurities allows for a number of simplifications that 
are not available for more general scattering mechanisms. 
First, the self-energy and the transport vertex are inde- 
pendent of momentum which allows us to employ shifts 
in the complex plane with help of which one can con- 
veniently express self-consistent quantities by their non- 
self-consistent counterparts. Second, the lack of momen- 
tum dependence leads to the lack of space non-localities 
of the scattering integral, therefore all virial corrections 
are described by the collision delay. Third, the momen- 
tum independence reflects that there is only a single scat- 
tering channel of s-symmetry. In general, different scat- 
tering channels have different collision delays, in the case 
of neutral impurities there is only a single collision delay 
what simplifies appreciably all related formulas. Fourth, 
due to the time-independence of the impurity potential, 
there are no gradient corrections to the scattering rate. 
Fifth, this time-independence also simplifies the energy 
conservation in collisions. Sixth, due to the absence of 
the dynamics of impurities, the virial corrections appear 
only in the scattering in. Here we have selected this sim- 
ple scattering on neutral impurities to have free hands to 
focus on details of the method. 

We have not discussed here consequences and interpre- 
tation of the virial corrections as it has been already done 
within the intuitive approach in Paper I. Our aim here is 
to confirm validity of this intuitive approach. 



14 



ACKNOWLEDGMENTS 



n B + n w f r . 



(A6) 



This work was supported from the Grant Agency of 
the Czech Republic under contract Nr. 202960098, the 
BMBF (Germany) under contract Nr. 06R0745(0), and 
the EC Human Capital and Mobility Programme. 



APPENDIX A: VIRIAL VERSUS 
QUASIPARTICLE CORRECTIONS IN 
EQUILIBRIUM 

The quasiparticle and the virial corrections enter the 
BE in different ways. From this point of view, they rep- 
resent independent corrections that can be treated sep- 
arately. On the other hand, for scattering by resonant 
levels, both corrections are of the same magnitude as 
it is demonstrated in Fig. 5 of Paper I. Striking simi- 
larity of their magnitudes raises a question up to what 
extent these two corrections are independent. To answer 
this question we briefly discuss equilibrium where one can 
benefit from well developed theory of quasiparticle .and 
virial corrections based on Green's functions.L3~E3't3 

In equilibrium, the local density of electrons is given 
by the spectral function as 



du 
~2tt 



dk 



■a(uj, k). 



(Al) 



In the limit of small scattering rates, the spectral function 
will be substituted from ( |l6| ) with (j> — and lowest order 
approximation of the scattering rate, 7 = 700- 

In the limit of small scattering rates we can easily sep- 
arate the quasiparticle contribution to the local density 



^-polc 



dk 



(A2) 



where f{k) = fFD^k), and the background contribution 



/duj , . f dk 
J (^ Re 



7oo M 



(u>-e k + iO) 2 



(A3) 



Using (p8|), ripoio can be split into the free part ( p3| ) and 
wave-function renormalization reduction 



where 



^pole — ^frcc 1 TT-wfr; 



dk d<T 00 
—f(k) 



(2tt 



(A4) 



(A5) 



From decompositions of t he density, n = n po i c 



ra fr , 



, and (A4) the correlated density results 



as a sum of the off-pole and wave-function renormaliza- 
tion parts 



Using the Kramers-Kronig relation for the real part of 
the self-energy 



coo M = Rc 



dE 



1 



2tt E-uj + iO 
, one finds that the correlated density reads 



dw dk ftuwv, 7ooM 

™corr = / — ,„ N , \JFD\U) — j(fcj) Re 



(A7) 



2tt (2tt) 3 



(A8) 



where the terms weighted by /f£>(w) and f(k) results 
from n ff and n w fr, respectively. Apparently there is a 
partial compensation of these contributions to n C orr- 

The quantum-mechanical expression (A8) c an b e also 
given the form of the semi-classical formula (102). To 
this end we reorganize n s starting from (A3), 



/du , . f dk 

-/™( W )7ooHRey 7^ 

= _ / £ /FD(w)7oo(w) c^ Re / 

du , , , ■. dgoo 
— / FD (w)7oo(w)-^— . 



n-off 



1 



) 3 (w-e fe + z0) 2 
dk 1 
(2tt) 3 Lu-e k +i0 

(A9) 



In the last line we have used that the integral over mo- 
mentum above defines a local element of Green's func- 
tion. 

To evaluate 700, we use 



700H = c(-2)Imi£,(u;), 



(A10) 



and the optical theorem [that follows from non-self- 
consistent form of (fll|)] 



Imt 



00 



t 



R\2 



Im 3oo- 



loo^yoo 



(All) 



If we express the local density of states in terms of the 
momentum integration 

- 2Img*(u;) = j A^o, ~ e P ), (A12) 

the off-pole contribution can be expressed in terms of the 
quasiparticle distribution, 



n-off 



dp du; R 2 dg 00 

f FD (Lu)\t oa (uj)\ 2n6(uj-e p ) 



(2tt) 3 2tt 

Jp_f (v )\ t R( e )|2 ^00 
(27r)3 /WI*oo^JI QuJ 



-n g and Finally, we substitute ( A13] ) into (| 



dp ef „\ ( dooo ,. R , 2 dgoa 

C|*00l -g^ 



(2tt) 



r/(p) 



du> 



dui 
(A13) 



(A14) 
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Form ula ( kXA ) is identical to the semi-classical expres- 
sion ( 102 ). To prove this claim we employ the derived 
optical theorem. 

It is advantageous to start from (|l02| ). First, we reor- 
ganize the ratio of the collision delay given by (pjJ) to the 
lifetime given by (59) as ~ = 2cImioo — * c (^oo — too) as 



— = ic(t 

T 



R 



n- 



1 c< 



i at* 



too 5w 



*oo duj 



2du [m+oo) 2 



Oct, 



III) 



duo 



l 00 UL 00 

too duj 



b 00 UL 00 

too duj 

L 00 UL 00 

t£ duj 



±R a.A 
L 00 ui 00 

too duj 



(A15) 



^From a non-self-consistent form of (El) one finds 



oz oo 



.2dg 



no 



duj 



(A16) 



which substituted into ( |A15| ) provides the derived optical 
theorem 



Oct, 



mi 



duj 



c\t 



R |2 



on 



du 



(A17) 



The l.h.s is exactly t he b racket in ( A14 ). Thus the 
semi-classical formul a (|102| ) is equivalent to the quantum- 
mechanical formula ( |A14| ). 

To summarize this appendix we would like to remind 
the most important point. Within Green's function, 
virial and quasiparticle corrections enter the density to- 
gether in an unresolved form, but the optical theorem 
and the derived optical theorem can be used to separate 
them and express virial corrections in terms of the colli- 
sion delay. 
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